Bistability and chaos in Taylor-Green dynamo 
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Using direct numerical simulations we study dynamo action under the Taylor-Green forcing with 
Prandtl number less than one. We observe bistability with a weak magnetic field branch and a strong 
magnetic field branch. Both the dynamo branches undergo subcritical dynamo transition. We also 
observe host of dynamo states including constant, periodic, quasiperiodic, and chaotic magnetic 
fields. One of the chaotic state originates through a quasiperiodic route with phase locking, while 
another chaotic attractor appears to follow Newhouse-Ruelle-Takens route to chaos. We also observe 
intermittent transitions among quasiperiodic and chaotic states for a given Taylor-Green forcing. 

PACS numbers: 91.25.Cw, 52.65.Kj, 47.20.Ky 



I. INTRODUCTION 



Dynamo theory has been applied to explain the gener- 
ation and properties of the magnetic fields present in ce- 
lestial bodies . In this mechanism, a small magnetic 
field fluctuation is amplified by the currents induced by 
the motion of the conducting fluid. One of the important 
problems of dynamo research is the nature of "dynamo 
transition" from pure fluid state to dynamo state. The 
linearized magnetohydrodynamic (MHD) equations yield 
zero magnetic held in the steady state. Therefore, the 
dynamo transition is through a nonlinear instability. As 
a result, the nature of transition is very rich exhibiting 
host of interesting behaviour including subcritical bifur- 
cations d, Q, variety of dynamo states d, 0], multiple 
coexisiting attractors Q etc. It has been found that 
the nature of dynamo onset depends critically on vari- 
ous system properties, e.g., the magnetic Prandtl num- 
ber (ratio of the kinematic viscosity and the magnetic 
diffusivity), forcing function, system geometry, rotation 
frequency, etc. Note that the magnetic Prandtl num- 
ber {Pm) of liquid metals, and of the convective fluids of 
Earth's outer core and the Sun is very small (of the order 
of 10~^), while that of the intergalactic medium is very 
large (around 10^*). In the present paper we explore the 
nature of dynamo transition for a Prandtl number less 
than unity, with a hope that it could be a representative 
of low-Prandtl number dynamo. 

One of the strong motivations of the dynamo research 
is to understand the nature of geodynamo and solar dy- 
namo. To this aim, scientists attempt to address the 
possible dynamo states, and the extent of the induced 
magnetic field in these systems. The magnetic Reynolds 
number Rm (ratio of magnetic advection and magnetic 
diffusion) of the geodynamo is estimated to be around 
125 which is somewhat near the dynamo transition, 
so we expect our present work on dynamo transition to 
be relevant for geodynamo. The ratio of the magnetic en- 
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ergy and the kinetic energy for the geodynamo is large [sj , 
hence it is called strong field dynamo. Roberts (Qj] argued 
that rotating sphere with magnetoconvection could show 
bistability with a weak field branch and a strong field 
branch. Kuang and Bloxham [lo| numerically simulated 
convective dynamo in spherical geometry with rotation 
and observed weak-field dynamo solution in a simplified 
system and strong-field dynamo solutions in a more re- 
alistic system. In the present paper we will show that a 
similar bistability is exhibited in the box geometry with- 
out rotation or convection. We also observe that the 
transition for both the branches are subcritical. 

The magnetic field of the Earth is temporally and spa- 
tially chaotic Q. The solar magnetic field is spatially 
random, but the sun-spot cycle appears to indicate that 
the primary dipolar field is quasiperiodic [llj . In the von 
Karman sodium (VKS) experiment, Monchaux et al. @ 
and Petrelis at al. reported various dynamo states 
including constant, periodic, quasiperiodic, and chaotic 
magnetic fields. In the present paper we simulate dynamo 
transition for a box geometry, and report various dynamo 
states (including chaos) similar to those obtained in the 
VKS experiment. 

We focus on the behaviour of Taylor-Green (TG) dy- 
namo for Pm — 0.5 fluid in order to probe dynamo 
behaviour for low magnetic Prandtl number (low-Pm) 
regime. Numerical simulation of very low-Pm fluid 
(around 10"'^ corresponding to liquid metals) is very dif- 
flcult since the corresponding Reynolds number for dy- 
namo transition for such fluid is more than 10^ 6J; simu- 
lation of such high Reynolds number would require grid 
resolution far more than presently achievable in the mod- 
ern supercomputers. We adopt TG forcing since it has 
somewhat similar flow structure as the VKS experiment. 

Earlier, Nore et al. (Tsl si mulated TG dynamo for Pm 
near unity. Ponty et al. 1J| studied dynamo mechanism 
for low magnetic Prandtl numbers (Pm « 1 to 10~^) 
using hyperdissipative parameters and large-eddy simu- 
lations, and observed a sharp increase in the critical mag- 
netic Reynolds number (threshold of dynamo transition) 
with the decrease of Pm. Ponty et al. observed a sub- 
critical dynamo bifurcation for the TG dynamo in the 
low-Pm regime. Mininni et al. (isj and Yadav et al. Q 
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studied the energy transfers, and the geometry of the ve- 
locity and the magnetic field structures. Yadav et al. p/j 
also observed a supercritical pitchfork bifurcation for the 
dynamo transition for Pm = 1. They reported a large 
number of dynamo states including constant, periodic, 
quasiperiodic, and chaotic magnetic fields. Dubrulle et 
al. [16] investigated various bifurcations in both hydro- 
dynamic (with no magnetic field) and magnetohydrody- 
namic systems under TG forcing. Scientists have also 
studied dynamo behaviour for various kinds of forcing, 
e.g., Roberts, ABC, Ponomarenko, and random, as well 
as for different geometries, e.g., box, cylinder, sphere 
etc. Dynamo transition and subsequent dynamo 

states have also been studied using low-dimensional mod- 
els of dynamo fis'l. 

The outline of the paper is as follows: The numerical 
procedure is described briefiy in Section |TI1 Bifurcation 
analysis is presented in Section IIIII We analyze several 
chaotic windows and the routes to chaos in Section ITVl 
We conclude in Section IVl 



II. SIMULATION METHODOLOGY 

The governing equations of dynamo are same as those 
of MHD, which are 

atU+(u-V)u = -Vp + (J X B) + z/V^u + F, (1) 

dtB = V X (u X B) + r^V^B, (2) 

V-u = 0, (3) 

V-B = 0, (4) 

where u is the fiuid velocity, J is the current density, B 
is the magnetic field, p is the hydrodynamic pressure, v 
is the kinematic viscosity, rj is the magnetic diffusivity, 
and F is the external force field. The density of the fluid 
is assumed to be unity. The three important parameters 
related to dynamo instability are the magnetic Prandtl 
number Pm = v/rj, the Reynolds number Re = UL/v, 
and the magnetic Reynolds number Rm = UL/rj, where 
U and L are the large velocity scale and the large length 
scale, respectively. Note that Rm = Re x Pm, hence 
only two of the above three parameters are independent. 
In our study, we fixed the magnetic Prandtl number to 
0.5 with v = 0.1 and r] = 0.2. 

We solve the MHD equations (Eqs. (HJl])) numerically 
for a box geometry of size (2tt)^ with periodic boundary 
conditions in all the directions. We use a pseudospectral 
code TARANG [l^ to carry out our simulations. We ap- 
ply fourth-order Runge-Kutta scheme for time advance- 
ment with dynamically adjusted dt chosen using the CFL 
condition: dt = Ax/\/20E^ where Ax is the grid spacing 
and i?" is the total kinetic energy. On the velocity field 
we apply Taylor Green forcing 
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where Fq is the amplitude of the forcing and ko is the 
wavenumber. We set ko equal to 2. Note that the TG 
forcing has components only along x and y directions. 
Numerical simulations reveal that the TG forcing induces 
counter rotating eddies, and it mimics the flow structure 
of the VKS experiment qualitatively [6j . 

The box is discretizcd uniformly in all the directions, 
with most of the simulations on 64'^ grid. This resolution 
was sufficient to resolve our simulation near the dynamo 
transition, as demonstrated by the fact that the product 
of the Kolmogorov length and the largest wavenumber 
lies between 1.3 — 9 for all our runs. We also verified 
the grid-independence by performing few runs of a given 
dynamo state on 128^ and 64^ grids. All our simulation 
were dealiased using the 2/3 rule. We performed approx- 
imately 150 simulations for various forcing parameters 
(Fo — 1 : 46), and studied the global kinetic and mag- 
netic energies, as well as the amplitudes of the velocity 
and magnetic Fourier modes. The magnetic Reynolds 
number in our simulations ranged from around 3 to 90. 

The importance of large scale modes have been am- 
ply highlighted in dynamo literature. In our present 
dynamo simulation, as well as for Pm — 1 reported 
earlier by Yadav et al. Q, some of the magnetic and 
kinetic Fourier modes play a dominant role. For the 
magnetic Reynolds numbers employed in these simula- 
tions, the dominant velocity modes are (±2, ±2, ±2), 
(±4, ±4, ±4), (±4, ±4,0), and the dominant magnetic 
modes are (0,0, ±1), (0,0, ±2), (0,0, ±3), (±2,±2,=f3), 
(=F2, =f2, ±1). Here, the three arguments refer to x, y, 
and z components of a wavenumber in the Fourier space. 
The (±2, ±2, ±2) velocity mode contains maximum ki- 
netic energy due to the TG forcing at A: = (2,2,2). The 
above set of modes carry more than 95% of the total en- 
ergy. Note, however, that higher wavenumber velocity 
and magnetic Fourier modes become important in the 
turbulent dynamo. We also point out that the above 
modes play similar roles as the large-scale modes used in 
the earlier dynamo literature [H [13] ■ 

In the present paper we analyze the time series and 
phase space plots of the Fourier modes. This exercise pro- 
vides us with information about various dynamo states 
as well as the nature of bifurcations. Primary and sec- 
ondary instabilities can be conveniently illustrated using 
bifurcation diagrams. In the next section we will de- 
scribe bifurcation diagrams for Pm = 0.5 constructed 
using some of the dominant magnetic Fourier modes. 



III. BIFURCATION ANALYSIS & STRONG 
AND WEAK DYNAMO BRANCHES 

We construct a bifurcation diagram by using the time 
averaged values of the magnitude of the magnetic Fourier 
mode B(0,0, 1) for different values of Fq (see Fig. [T]). 
These values are computed after the system reaches a 
steady state. Since the Fourier amplitudes are generally 
complex, it is convenient to use the absolute value of a 
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FIG. 1: Bifurcation diagram: Time averaged values of 
|B(0, 0,1)1 for various forcing amplitude Fq. The figure il- 
lustrates various dynamo states with ND — no-dynamo state, 
FP = fixed point (constant in time), P — periodic state, QP 
= quasiperiodic state, and C = chaotic state. The inset for 
Fo = 14:16 illustrates a sudden jump and a hysteresis loop 
indicating that the transition is subcritical. Also, oval 'A' 
exhibits quasiperiodic route to chaos (see Subsection IIV A| l . 
and the rectangular box 'B' exhibits Newhouse-Ruelle-Takens 
route to chaos (see Subsection IIV Bl) . 



Fourier mode to depict the dynamo states. 

The inset of Fig. [1] shows a zoomed view of the dy- 
namo transition region near Fq = 15.8. We observe fluid 
or no-dynamo state for Fq < 15.8. At Fq = 15.8, we ob- 
tain a dynamo state with |B(0, 0, 1)| as well as the total 
magnetic energy {E^) showing a finite jump. However, 
when we use the dynamo state at Fq = 15.8 as our initial 
condition and gradually decrease the forcing amplitude, 
as shown in the inset of Fig. [1] the dynamo state con- 
tinues till Fq w 15.2, at which point there is a sudden 
jump to the fluid state. This feature of hysteresis for 
Fq = 15.2 : 15.8 demonstrates the subcritical nature of 
the transition, with a "subcritical pitchfork bifurcation" 
at X and a "saddle-node bifurcation" at Y (see the inset 
of the figure). These points are joined together using a 
dashed curve to depict the unstable branch. In a recent 
study, Krstulovic et al. [2l| also observed subcritical dy- 
namo transition in their simulation and low-dimensional 
model. In a planetary dynamo context, Kuang et al. [2^ 
argued that the sudden termination of Martian dynamo 
could be due to subcriticality. 

Another noteworthy feature during the transition at 
Fq — 15.8 is that the Reynolds number is sufficiently 
high (w 60) for the fluid to be temporally chaotic. How- 
ever, in the corresponding dynamo states at the same 
Fq, both the velocity and magnetic fields become con- 
stant in time, as shown in Fig. [5J This is due to the 
presence of the newly-born finite magnetic field. Morin 
et al. [2^ reported similar effects of magnetic field for 
spherical shell dynamo. After Fq 16.2, the dynamo 
solution bifurcates to periodic states, and subsequently 




F„ = 15.8 



500 



1500 
Time 



2500 



3500 



4500 



FIG. 2: The total kinetic energy (top panel) and the total 
magnetic energy (bottom panel) for a forcing amplitude of 
15.8 indicating a constant value for the steady state. 



to quasiperiodic and chaotic states as shown in Fig. [TJ 

For some of our dynamo runs, (|B(0, 0, 1) | is negligible, 
as shown by the dots on the x axis of Fig. [T] For these 
simulations, the magnetic mode B(0, 0, 2) becomes dom- 
inant; we denote the set of these states as the "B(0, 0, 2)- 
branch" . The other collection of states with B(0, 0, 1) as 
the most dominant mode is referred to as the "B(0, 0, 1)- 
branch" . We observe that the properties of these two 
branches are quite different. We illustrate this feature 
using another bifurcation diagram whose vertical axis 
is E^/E'^, the ratio of the magnetic energy and the ki- 
netic energy (shown in Fig. [S]). This bifurcation diagram 
reveals two distinct branches: the upper curve is the 
B(0, 0, 2)-branch, for which the magnetic field remains 
constant in time, and the lower one is the B(0,0, 1)- 
branch. The B(0, 0, 2) branch was constructed using the 
fixed point solution at _fo — 39.7, which is marked with 
small arrows in Figs. |3] and [1] as initial condition. Also, 
the ratio E''/E" for the B(0, 0, 2)-branch is larger than 
the corresponding ratio for the B(0, 0, l)-branch. Thus 
we address the former as a "strong-field branch", and 
the latter as a "weak-field branch". The above feature 
regarding the strength of the magnetic field is corrobo- 
rated in Fig. |4] where we plot the magnetic energy {E'') 
vs. the magnetic Reynolds number (Rm). The two dis- 
tinct branches described above demonstrate that the TG 
dynamo at Pm = 0.5 exhibits "bistability" , similar to 
those reported by Simitev and Busse [13] . 

The dynamo transition for the B(0, 0, 2)-branch differs 
significantly from that of the B(0, 0, l)-branch described 
earlier. For the B(0, 0, l)-branch, the jumps in the kinetic 
and magnetic energy shown in Fig.[5ja) is consistent with 
the subcritical nature of the transition described earlier. 
However, the B(0, 0, 2)-branch exhibits chaos at the dy- 
namo transition itself. The origin of the chaotic state 
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FIG. 3: Bifurcation diagram: Ratio of the total magnetic 
energy and tiie total kinetic energy {E^'/E") vs. the forcing 
amplitude (-fo)- Description of different symbols is provided 
in Fig. [T] The B(0, 0, 2)-branch was constructed using the 
dynamo state marked with an arrow as initial condition. 




FIG. 4: Bifurcation diagram: Total magnetic energy (E^) vs. 
the magnetic Reynolds number (Rm). Description of different 
symbols is provided in Fig. [1] 



becomes apparent when we decrease Fq from the fixed 
point solution. As shown in Fig. [5][b), the fixed point 
bifurcates to quasiperiodic (purple triangles of the fig- 
ure) and subsequently to chaotic state (red squares) as 
the Fq is decreased. The route to chaos for the B(0, 0, 2)- 
branch may be similar to those observed by Pal et al. [25j 
for zero-Prandtl number convection. The quasiperiodic 
and the chaotic states are clustered in Fig. [3] near the 
transition. 

Earlier, the weak field branch and strong field branch of 
dynamo have been reported in the context of convection- 
driven hydromagnetic dynamos [1, [IS IHj 113 • Pierre [ll] 
and Stellmach and Hansen [2^ numerically demonstrated 
this phenomenon in Childress-Soward dynamo where 
they reported that the Strong-field branch goes even 
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FIG. 5: A plot of the total kinetic energy (-B") and the total 
magnetic energy {E'') vs. the Fo near the dynamo transition. 
Description of different symbols is provided in Fig. [T] (a): 
The B(0, 0, l)-branch portrays the subcritical nature of the 
transition, (b): The B(0, 0, 2)-branch illustrates transition 
from a fixed point state to quasiperiodic (purple triangles) and 
subsequently to chaotic state (red square) as Fo is decreased. 



further back than the weak-field branch in the parame- 
ter space, analogous to the B(0, 0, 2)-branch reported in 
this paper. Similar behaviour has also been reported by 
Sreenivasan and Jones j30| in a rapidly rotating spher- 
ical shell dynamo. It is interesting that our system in 
the box geometry, without any convection or rotation, 
exhibits weak and strong dynamo branches. 

In the next section we will focus on several chaotic 
dynamo states and study routes to chaos for these states. 



IV. ROUTES TO CHAOS 

The phase space portraits drawn using the Fourier 
modes reveal multiple windows of chaos near the dynamo 
transition itself. In the present paper we describe routes 
to chaos for two chaotic windows: the oval-shaped en- 
closure 'A' and the rectangular-shaped box 'B' shown 
in Fig. [T] In the following discussion we show that the 
origin of chaos for the dynamo states of box 'A' follow 
a quasiperiodic route to chaos through phase locking, 
while the corresponding route to chaos for the box 'B' 
is through the Newhouse-Ruelle-Takens scenario. 

We study the routes to chaos using some of the dom- 
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FIG. 6: Phase space projections on the |B(0, 0, 1)|-|B(0, 0, 3)| 
plane demonstrating a quasiperiodic route to chaos through 
phase-locking. The forcing amplitudes (Fo) for different at- 
tractors are marked with arrows. _Fo — 41, 42, 43 (pink), 
Fo = 43.85 (green), and Fo = 44.75, 46 (red) correspond to 
quasiperiodic, phase-locked, and chaotic states, respectively. 



inant modes of magnetic field, viz. B(0, 0, 1), B(0, 0, 2), 
and B(0,0, 3). We use time series, phase-space projec- 
tions, and Poincare sections for this study |31| . 
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FIG. 7: Poincare sections for some of the attractors of Fig. [6] 
with |B(0, 0,3)1 = 1.15 as the Poincare plane, (a): quasiperi- 
odic attractor for Fq = 43, (b): phase- locked attractor for 
Fo = 43.85, (c,d): chaotic attractors for Fo = 44.75 and 46. 



B. Newhouse-Ruelle-Takens scenario 



Quasiperiodic route to chaos through 
phase-locking 



In Fig. |6] we illustrate the phase space projections of 
dynamo states on |B(0, 0, 1)|-|B(0, 0, 3)| plane for the 
forcing range of Fq — 41:46 (corresponding to the oval- 
shaped enclosure 'A' of Fig. [T]). For Fq = 41, 42, and 43, 
the system is quasiperiodic since the phase space projec- 
tion is densely filled up. The approximate values of the 
two incommensurate frequencies for the Fq = 43 dynamo 
state are 0.0165 and 0.0208 (for |B(0,0,1)| time series). 
For Fq — 43.85, the system becomes periodic or phase- 
locked. The time period of the periodic orbit is relatively 
large. The emergence of periodic orbit after quasiperiodic 
solutions is called "phase locking" [3l| . A subsequent in- 
crease of Fo leads to a chaotic state, as evident from 
the phase space projections for Fq = 44.75 and 46. The 
nature of attractors is corroborated by the Poincare sec- 
tions for Fo 43, 43.85, 44.75, and 46 presented in Fig. [7] 
(a, b, c, d), respectively. These Poincare sections were 
obtained by using the |B(0,0, 3)| — 1.15 as the Poincare 
intersection plane for the phase space trajectories in the 
subspace of |B(0, 0, 1)|-|B(0, 0, 2)|-|B(0, 0, 3)|. Thus the 
route to chaos for the chaotic trajectories in the oval-box 
'A' of Fig.[T]is through the quasiperiodic route via phase 
locking, first reported for circle map jSlj]. 



We observe an interesting set of dynamo states in 
the rectangular-box 'B' of Fig. [TJ Fig. IUJa) illustrates 
the phase space projections on the |B(0, 0, 1)|-|B(0, 0, 3)| 
plane for Fq = 39.3 and 39.464, corresponding to the 
two purple triangles in the rectangular-box 'B' of Fig. [T] 
Also, the power spectral density plot (PSD) of the 
|B(0, 0,1)1 time series of these two dynamo states is 
shown in Fig. ini^a,b). The PSD reveals that the state 
at Fo = 39.3 contains two incommensurate frequencies 
/i and /2 (Fig. [IJa)), while the state at Fo = 39.464 
has three incommensurate frequencies /i, /2, and /a 
(Fig. Eljb)). Thus the corresponding dynamo states re- 
side on "2-torus" (T^) and "3-torus" (T^) respectively in 
the subspace. We expect that a further increase of Fo 
should push the system to chaos following a Newhouse- 
Ruelle-Takens scenario. However, we have not yet found 
the corresponding chaotic attractor. Recently, Stefani 
et al. [32] have observed similar quasiperiodic routes to 
chaos in a dynamo model. 

Instead of transition from to a chaotic dynamo 
state, a very small increase of Fq pushes the system to 
a new chaotic attractor, whose span is much larger than 
those of the or T^. In Fig.|8l[b), we illustrate the larger 
chaotic attractor obtained for Fq = 39.51, corresponding 
to the red square in the rectangular-box 'B' of Fig. [T] 
In Fig.mjb), or states reside in the circled region. 
The large attractor for Fo = 39.51 has a different ori- 
gin as the PSD shown in Fig. ^c) is quite different from 
those in Fig. [5] (a, b). The origin of the larger attractor 
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FIG. 8: Phase space projections on the (|B(0, 0, 1)|- 
|B(0, 0,3)1) plane. Panel (a) contains a stable 2-torus and 
a 3-torus quasiperiodic state and panel (b) shows a close ly- 
ing (in parameter space) large chaotic attractor. The corre- 
sponding forcing amplitude is marked with arrow. The circle 
drawn in (b) portrays a rough phase space span of the attrac- 
tors shown in (a). 



is not well understood, however, it is possibly through a 
"crisis" . Note that a large number of attractors exists for 
this range of Fq, as evident from the Fig. [TJ Chaos can 
emerge due to the intersections of multiple attractors or 
their basins of attraction [3l|. A detailed investigation 
of these issues is beyond the scope of this paper. 

We performed dynamo simulations near Fq = 39.51 
and observe interesting features involving intermittent 
transitions between three attractors; this dynamics will 
be described in the next subsection. 
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FIG. 9: Power Spectral Density (PSD) plot of the dynamo 
states shown in Fig. [S] The dominant frequencies are depicted 
in the figures with arrows, (a): For Fq = 39.3, the state has 
two incommensurate dominant frequencies; (b): For Fq = 
39.464, the state has three incommensurate frequencies; (c): 
A dense power spectrum corresponding to the large chaotic 
attractor portrayed in Fig.igb). 



CONCLUSIONS AND SPECULATIONS 



C. Intermittent transitions between various 
attractors 

For the forcing amplitude Fq — 39.4658 we observe 
intcrmittcncy. A long time series for Fq — 39.4658 
simulation is illustrated in Fig. llOf a). which illustrates 
that the system makes intermittent transitions among 
three attractors, '2F', '3F', 'C, shown below the time 
series. Note that the 2F, 3F, and C attractors are qual- 
itatively similar to the state shown in Fig. [5] obtained 
for Fq = 39.3,39.464, and 39.51, respectively. Thus, at 
Fq ~ 39.4658 the system appears to hop over various 
attractors. The intermittent transition among various 
attractors described above is similar to "intermittency" 
in which a system switches between an "ordered" state 
and a chaotic state for a single parameter. In Fig.fTOlthe 
quasiperiodic states are the ordered state. The system 
has been evolved for 50000 eddy turnover times, and only 
a part of this time series has been shown in Fig. rTOT a). 
Hence the observed phenomena is not a transient. The 
large fluctuations for time near 21250 in Fig. [TUf a) are 
transient fluctuations that settle down quickly. 



Our Taylor-Green dynamo simulations for Prii = 0.5 
near the dynamo transition reveal interesting proper- 
ties. We observe bistability, with a weak magnetic field 
branch and a strong magnetic field branch. The domi- 
nant Fourier modes for the two branches are B(0,0, 1) 
and B(0,0, 2) respectively. Both of these branches have 
subcritical origin. Qualitatively, these branches resemble 
the weak- and strong-field dynamo actions in rotating 
magnetoconvection reported earlier [l^] . Note that 
the subcritical dynamo transition for Pm = 0.5 is in 
sharp contrast with the TG dynamo for Pm = 1 which 
shows a supercritical pitchfork bifurcation at the transi- 
tion 0. This observation is consistent with the recent 
work by Krstulovic et al. [2l[ . 

Our simulations also reveal various kinds of dynamo 
states including constant, periodic, quasiperiodic, and 
chaotic magnetic fields. We analyzed two chaotic win- 
dows among several ones observed in our simulations. 
Chaos for these windows arise through quasiperiodic 
route. One of them undergoes a phase-locking, while the 
other one appears to be through the Newhouse-Ruelle- 
Takens route. We also observe an intermittent transi- 
tions between quasiperiodic {T^,T^) and chaotic attrac- 
tors. Various dynamo states reported in our simulations 
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FIG. 10: Time evolution of a dynamo state at _Fo ~ 39.4658 showing intermittent transitions among three attractors: quasiperi- 
odic state T^(b), quasiperiodic state T^{c), and chaotic state(d). The subfigures (b,c,d) are the phase space projections at 
three different time intervals during the evolution. 



are quite similar to those observed in the VKS experi- 
ment [1]. 

We performed our simulations on a periodic box, which 
is an idealized geometry. Yet, we observe similarities with 
the VKS experiment and geodynamo simulations. This 
could be due to similarities in the inherent nonlinearities 
in the system (independent of the geometry etc.). Also, 
the behaviour of dynamo transition and dynamo states 
for Pm = 0.5 discussed in the present paper differs sig- 
nificantly from that for Pm = 1 studied by Yadav et 
al. Extensions of the present work to more realis- 
tic geometry and lower magnetic Prandtl numbers would 



yield interesting insights into this challenging field. 
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